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ABSTRACT 

We present a detection of [CI] 3 Pi— 3 P emission in the lensed quasar 
APM 08279+5255 at z=3.91 using the IRAM Plateau de Bure interferometer. 
The [CI] line velocity and width are similar to the values of previously detected 
high-J CO and HCN lines in this source, suggesting that the emission from all 
of these species arises from the same region. The apparent luminosity of the 
[CI] line is L' cl = (3.1 ± 0.4) x 10 10 K km s^ 1 pc 2 , which implies a neutral 
carbon mass, Mci = (4.4 ± 0.6)m" 1 x 10 7 M , where m is the lensing magni- 
fication factor. The [CI] line luminosity is consistent with the large molecular 
gas mass inferred from the nuclear CO line luminosity (~10 11 m~ 1 M ). We also 
present an upper limit on the H2O lio— loi nne luminosity in APM 08279+5255 
of, L^q < 1.8 x 10 10 K km s" 1 pc 2 (3-ct). 

Subject headings: - galaxies: active - galaxies: high-redshift - quasars: emission 
lines - quasars: individual (APM 08279+5255) - galaxies: ISM 



1 Based on observations carried out with the IRAM Plateau de Bure Interferometer. IRAM is supported 
by INSU/CNRS (France), MPG (Germany) and IGN (Spain). 
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1. Introduction 

The ultraluminous quasar APM 08279+5255 at z=3.91 (Irwin et al. 1998) is strongly 
gravitationally lensed (Ledoux et al. 1998; Ibata et al. 1999; Egami et al. 2000) and appears 
to be one of the most luminous objects in the Universe, with an apparent infrared luminosity 
of ~10 15 L (Lewis et al. 1998). Observations over a wide range of wavelengths (e.g., Downes 
et al. 1999; Ellison et al. 1999; Gallagher et al. 2002; Soifer et al. 2004; Wagg et al. 2005) have 
shown that APM08279+5255 contains an active nucleus and likely a starburst component. 

The combination of extreme intrinsic luminosity and amplification by strong gravita- 
tional lensing has allowed APM 08279+5255 to be detected in many rotational lines of CO 
(J=4-3 and J=9-8 Downes et al. 1999, hereafter D99; J=l-0 and J=2-l Papadopoulos et al. 
2001; Lewis et al. 2002 ), including the unique detection in the high excitation J=ll-10 line 
(A. Weiss et al., in preparation). From high angular resolution observations, D99 suggest 
that the high-J CO emission arises in a warm, dense circumnuclear disk of sub-kiloparsec 
size. The detection of both HCN J=5-4 and HCO + J=5-4 line emission support the idea 
that the molecular gas is dense (Wagg et al. 2005; Garcia-Burillo et al. 2006) 

Given the wealth of molecular lines detected in APM 08279+5255, observations of ad- 
ditional diagnostic lines can provide further constraints on the physical conditions in this 
high-redshift object. One important tracer of the dense neutral gas within a galaxy's in- 
terstellar medium is atomic carbon, in particular the [CI] 3 Pi— 3 Po line with rest frequency 
492.161 GHz (Gerin & Phillips 1998, 2000). This line has not been observed widely in local 
galaxies due to poor atmospheric transmission near this frequency; however, this line is red- 
shifted into better atmospheric windows for objects at redshifts, z > 1. The [CI] 3 Pi— 3 P 
line provides a complimentary probe of the physical conditions of dense, neutral gas, and the 
[CI] /CO ratio may also provide information on the role of X-rays in the molecular excitation 
(Maloney et al. 1996). In addition, Papadopoulos et al. (2004a, b) argue that this line can 
provide a measure of the total molecular hydrogen (H 2 ) gas mass, independent of CO line 
emission. 

To date, [CI] emission has been detected in four objects at z > 2, most of which are 
thought to be gravitationally lensed. The 3 Pi— 3 Po line was first detected in H1413+117 
at z=2.6 (the "Cloverleaf " ; Barvainis et al. 1997), and more recently the higher excitation 
3 P2— 3 Pi line was also detected in this object (Weiss et al. 2003). Other high-redshift objects 
with [CI] 3 Pi— 3 Po line detections include IRAS F10214 at z=2.3 (and an upper-limit on the 
[CI] 3 P 2 - 3 Pi line intensity is reported by Papadopoulos 2005), SMM J14011+0252 at z=2.6 
(Weiss et al. 2005), and PSS 2322+1944 at z=4.1 (Pety et al. 2005). 

Another important species that acts as a major coolant of dense gas is water. In general, 
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H 2 lines at local velocites are not accessible from the ground, but they have been observed 
from space and are ubiquitous in dense molecular cloud cores in the Galaxy (Ashby et al. 
2000; Snell et al. 2000). At high-redshift, a tentative detection of the H 2 2 n — 2 02 line, with 
rest frequency 752.033 GHz, was reported in IRAS F10214 at z=2.3 (Encrenaz et al. 1993; 
Casoli et al. 1994). At more modest redshifts, the fundamental transition of ortho- water, 
H 2 lio — loi has been detected in absorption towards B0218+357 at z=0.685 (Combes & 
Wiklind 1997). For APM 08279+5255, the lower excitation H 2 li -l i line with rest 
frequency 556.936 GHz is redshifted into the 3 millimeter atmospheric window. 

Here we report a detection of [CI] 3 Pi— 3 Po line emission in APM 08279+5255 at z=3.91, 
and an upper-limit on H 2 lio— loi line emission, obtained with the IRAM Plateau de 
Bure Interferometer (PdBI). Throughout this paper we adopt a A-dominated cosmology: 
H = 70 km s" 1 Mpc-\ fi A = 0.7, and tt m = 0.3 (Spergel et al. 2006). 

2. Observations 

We used the IRAM PdBI to search for the [CI] 3 Pi- 3 P and H 2 li -l i lines in 
APM 08279+5255, redshifted to the 3-millimeter band, on six dates in 1999 and 2001. All of 
the observations were made in compact configurations of the four or five available antennas. 
The receivers were tuned to 100.216 GHz for the [CI] line, and 113.406 GHz for the H 2 line. 
Spectral correlators covered a velocity range of ~1500 km s -1 for [CI], and ~1300 km s -1 
for H 2 0. The phase center was offset by 1.25" from the CO peak position in D99. Baseline 
lengths ranged from 17 to 81 meters, and the synthesized beam sizes were 7'.'2 x 6'/0 (position 
angle 89°) for the [CI] image, and 5'.'9 x 4'.'4 (position angle 82°) for the H 2 image. The 
nearby quasar 0749+540 was used for complex gain calibration. The flux scale was set 
using standard sources including MWC349, CRL618, 3c345.3 and 0923+392, and should be 
accurate to better than 20%. The on-source integration time for the [CI] and H 2 images 
were equivalent to 7.3 hours and 3.1 hours, respectively, with the full six antenna array. The 
resulting [CI] image has an rms noise of 0.36 mJy beam -1 in a 150 km s -1 channel, and the 
H 2 image has a higher rms noise of 0.96 mJy beam -1 in a 130 km s -1 channel, reflecting 
lower atmospheric transmission near the edge of the 3-millimeter window. 

3. Results 

Figure 1 shows images of the [CI] line and continuum emission over the full velocity 
range observed. The line emission is clearly detected in several velocity bins and appears 
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spatially unresolved. The position of peak emission is consistent with that found previously 
for the dust continuum, CO lines, and HCN line. Figure 2 shows the spectra of the [CI] and 
H 2 lines at this position. There is no evidence for significant H 2 line emission. For the 
[CI] spectrum, we estimate the continuum level at 100.2 GHz using the "line-free" channels 
to be S'iooghz = 1.18 ± 0.18 mJy. For the H2O spectrum, we estimate the continuum level 
at 113.4 GHz using all of the channels to be Susqkz = 2.43 ± 0.48 mJy. These values 
are consistent with the 93.9 GHz continuum flux of 1.20 ± 0.30 mJy measured by D99 and 
a thermal spectrum. A Gaussian fit to the continuum subtracted [CI] 3 Pi— 3 P spectrum 
yields a peak of 2.20 ±0.51 mJy, central velocity v = 117 ±28 km s" 1 (z=3.9130 ± 0.0005), 
and AVfwhm = 386 ± 67 km s~ x . The integrated intensity of the [CI] 3 Pi— 3 P line is 
0.93 ±0.13 Jy km s _1 , which implies a line luminosity L' cl = (3.1 ±0.4) x 10 10 K km s" 1 pc 2 , 
following Solomon et al. (1992). Table 1 lists the fitted and derived [CI] 3 Pi— 3 P line 
parameters. 

We place an upper limit on the H2O lio— loi line emission from the noise around 
the mean in the spectrum, assuming AVfwhm = 450 km s _1 like the CO, HCN, and [CI] 
lines. The 3-er upper limit to the integrated intensity is 3 • (AVFWHM/'AV'chan) 1 ^ 2 • ^chan — 
0.70 Jy km s _1 , which implies an upper limit to the H 2 line luminosity -^h 2 o < 1-8 x 
10 10 K km s' 1 pc 2 . 

4. Discussion 

4.1. Emission Region and [CI]/CO Luminosity Ratio 

The fitted [CI] line center and width are compatible with previous observations of the 
HCN J=5-4 and HCO± J=5-4 lines (Wagg et al. 2005; Garcia-Burillo et al. 2006), though 
offset by ~ 120 km s _1 from the mean redshift determined from the high-J CO lines (D99; 
Weiss et al. 2006, in preparation). While the line center difference between the [CI] and CO 
is formally significant at the ~4<r level, the true significance is likely lower, since the line 
profiles from this complex source are almost certainly not described accurately by a single 
Gaussian. Given this uncertainty, the modest signal-to-noise ratios, and the lack of any 
resolved velocity structure in the observed line profiles, we simplify our analysis by assuming 
the [CI] emission region is cospatial with the CO line emission region, such that the conditions 
derived for a model with a single physical component apply to both species. We stress that 
even if all of the observed lines do arise from a single physical component, the fitted Gaussian 
centers and widths of lines with different excitation properties will not necessarily match each 
other perfectly, given the complexity of the underlying source structure and optical depth 
effects. Data of higher quality and from more transitions will be required to relax the single 
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component assumption and to justify more sophisticated modeling. 

The gas in the single component region is warm and dense, and is likely to be a cir- 
cumnuclear disk of sub-kiloparsec size. The CO emission clearly has been magnified by 
gravitational lensing, and models suggest a magnification factor in the range 3 to 20 (D99, 
Lewis et al. 2002), or even higher in the models by Egami et al. (2000). The exact mag- 
nification factor depends critically on the intrinsic size of the emitting region. If the [CI], 
CO, and HCN emission are truly co-spatial, however, then the magnification factors are the 
same, and line luminosity ratios will not suffer any biases due to differential lensing effects. 

The luminosity ratio between the [CI] and CO lines is a potential diagnostic of the 
gas excitation and carbon chemistry. To compare the luminosities of [CI] and CO in 
APM 08279+5255, we use the nuclear component of the CO J=l-0 emission (Lewis et al. 
2002), which gives Lf CI] /L' co(1 _ 0) = 0.23 ± 0.06. This value falls within the range 0.2 ± 0.2 
found for local galaxies (Gerin & Phillips 2000), and also the range ~ 0.15 - 0.32 for the 
other four high-redshift objects that have been detected in both [CI] and CO emission. This 
ratio does not appear to have any strong dependence on environment, either in the local 
galaxy sample, or at high-redshifts. Given the significant uncertainties in the [CI] /CO lu- 
minosity ratio for any individual object, it is not possible to place strong constraints on the 
contribution of X-ray irradiation to the molecular gas excitation. 

We can estimate the [CI/CO] abundance ratio from the observations, adopting the 
physical conditions derived from multi-transition CO observations (Wagg et al. 2005; A. 
Weiss et al. in preparation) and a simple radiative transfer model. Using the RADEX 2 LVG 
code (Schoier et al. 2005), we find that an abundance ratio, [CI/CO] ~ 0.6, reproduces the 
observed [CI] 3 Pi— 3 Po to CO J=l-0 line luminosity ratio. This ratio is similar to that found 
in the nearby starburst galaxy M82 ([CI/CO] ~ 0.5; Schilke et al. 1993, White et al. 1994). 



4.2. Neutral Carbon Mass 

The luminosity of the [CI] 3 Pi— 3 P line can be used to derive the neutral carbon mass 
(Weiss et al. 2005), as 

M [CI] = 5.706 x 10-*9iM e ^/T< xL > ci] [Mq]> (1) 

where Q{T ex ) ~ l + 3e~ Tl / T - + 5e~ T2 ^ (T\ = 23.6 K and T 2 = 62.5 K are the energies of the 
3 Pi and 3 P 2 levels above the ground state) is the partition function. The only high-redshift 



2 http://www. strw.leidcnuniv.nl/~moldata/radex. html 
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object for which both the [CI] 3 Pi— 3 P and [CI] 3 P 2 — 3 Pi lines have been measured is the 
"Cloverleaf" quasar, where T ex = 30 K (Weiss et al. 2003, 2005). Without a measurement of 
the higher [CI] 3 P2— 3 Pi line in APM 08279+5255, the value of T ex is uncertain. However, 
Weiss et al. (2005) show that for T ex > 20 K and LTE, the estimated neutral carbon mass is 
insensitive to the assumed excitation temperature. Adopting a nominal value of T ex = 80 K 
for the [CI] 3 Pi— 3 Po line in APM 08279+5255 (from the derived molecular gas and dust 
temperatures), the neutral carbon mass within the nuclear region of APM 08279+5255 is 
Ma = (4.4 ± Q.Q)m~ l x 10 7 M , where m is the lensing magnification factor. If T ex were 
20 K, then the estimated [CI] mass would be ~30% lower. 



4.3. Molecular Gas Mass 

The luminosity in the CO J=l-0 line is the traditional tracer of molecular gas (= H2) 
mass in the Galaxy and in a variety of extragalactic environments. For nearby ultraluminous 
infrared galaxies (ULIRGs; Sanders & Mirabel 1996; Sanders et al. 1988), which are believed 
to be analogues to the currently detectable luminous high-redshift galaxies, the empirical 
scaling factor between CO J=l-0 luminosity and H 2 gas mass is ~1 M Q (K km s^ 1 pc 2 ) -1 
(Downes & Solomon 1998). If we apply this factor to the nuclear CO J=l-0 line emission 
in APM 08279+5255 (Papadopoulos et al. 2001; Lewis et al. 2002), then the molecular gas 
mass implied is ~ 1.5m -1 x 10 11 M Q . 

Following Papadopoulos et al. (2004a), the small scatter in the observed [CI] /CO lu- 
minosity ratio suggests that [CI] 3 Pi— 3 P emission may provide an independent estimate 
of molecular gas mass. Papadopoulos et al. (2004a) relate the integrated [CI] line intensity, 
S[ci], to molecular (= H 2 ) gas mass, as 

M H2 _, CI] = 13 T5,^L (^f ( I ^)- 1 0rolp -| f _ [Ms] , (2) 

where X CI is the [CI]-to-H 2 abundance ratio (we assume X CI = 3 x 10~ 5 ), A w is the Einstein 
A-coefficient (A w = 7.93 x 10~ 8 s _1 ), and Q w is the excitation factor which depends on the 
kinetic temperature, and density of the gas (we adopt Q w = 0.5). The cosmology dependence 
is included through the luminosity distance, Di. Using the observed integrated intensity of 
the [CI] 3 Pi- 3 P line in APM 08279+5255 gives M H2 _ [C i] = (2.7 ± O^m" 1 x 10 11 M , in 
good agreement with the molecular gas mass estimated from the CO J=l-0 line. 

Table 2 presents the molecular gas mass estimates derived from [CI] and from CO for the 
five high-redshift objects where both lines are detected, adjusted to a common cosmology. 
The [CI] 3 Pi- 3 P line integrated intensities for IRAS F10214, SMM J1401 1+0252, and 
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H1413+117 are taken from Weiss et al. (2005), and for PSS 2322+1944 from Pety et al. 
(2005). The CO J=l-0 line has been detected only in APM 08279+5255 and PSS 2322+1944 
(Carilli et al. 2002). For the other objects, which are not detected in the CO J=l-0 line, 
we use observations of the CO J=3-2 line (IRAS F10214; Solomon et al. 1992, H1413+117; 
Weiss et al. 2003, SMM J14011+0252; Downes k Solomon 2003), which is approximately 
equal to the luminosity in the CO J=l-0 line at high-redshifts for warm, dense gas (Solomon 
et al. 1992). Though the uncertainties are large, there is very good agreement between the 
molecular gas mass estimates from the two tracers. 



4.4. H 2 

We did not detect any significant H2O lio — loi nne emission in APM 08279+5255. We 
can estimate the H2O line luminosity expected for the circumnuclear region for a nominal 
water abundance, again adopting the previously derived physical conditions for the region. 
In Galactic molecular clouds, a typical ortho-H 2 abundance is 4.5 x 10~ 9 relative to H 2 
(Snell et al. 2000; Ashby et al. 2000). For T kin = 80 K, n H2 ~ 40,000 cm" 3 , an LVG 
calculation gives N(R 2 0)/Av = 8.0 x 10 12 cm- 2 (km s -1 ) -1 , L H2 o = 1.1 x 10 7 K km s^ 1 pc 2 , 
more than two orders of magnitude below the upper limit. In strongly shocked regions, the 
H 2 abundance has been observed to be enhanced by more than an order of magnitude, 
but even if that were the case, calculations show that the H 2 lio— loi li ne emission would 
remain well below the achieved detection threshold. 



5. Summary 

We detected [CI] 3 Pi— 3 Po line emission in the ultraluminous quasar APM 08279+5255 
at z=3.91. The [CI] 3 Pi— 3 Po line width and center are similar to those of previously de- 
tected millimeter CO and HCN lines. Estimates of the molecular gas mass based on the 
[CI] 3 Pi— 3 Po line and based on CO lines yield similar results, and there is no evidence for 
any substantial contribution from a primarily atomic medium. Though the observational 
uncertainties remain large, this seems to be the case for the five high-redshift sources with 
reported detections of both CO and [CI] emission lines. The [CI] 3 Pi— 3 Po line may be 
a valuable alternative probe of molecular gas mass in systems lacking high excitation CO 
line emission, where the molecular medium is more diffuse and cooler (T^n < 20 K,and 
riH 2 i$ 10 4 cm -3 ; Papadopoulos et al. 2004a). The Atacama Large Millimeter Array will 
greatly expand the number of high-redshift sources accessible in [CI] and CO emission, and 
will have the capability to spatially resolve the emission to show directly the extent of the 
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gas. 
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Table 1. [CI] 3 Pi- 3 P line parameters for APM 08279+5255. 



[CI] 3 P r 


— 3 P peak: 


2.20 ±0.51 mJy 


[CI] 3 P r 


— 3 P AVf WH m : 


386 ± 67 km s" 1 


[CI] 3 P r 


- 3 Po v a : 


117 ±28 km s" 1 


S[ci]- 




0.93 ±0.13 Jy km s" 1 


Tl b . 

L [CI] ■ 




3.1 ±0.4 x 10 10 K km s" 1 pc 2 



a velocity with respect to z=3.911, determined from CO 
lines. 

b The [CI] 3 Pi— 3 P luminosity is not corrected for lensing 
magnification. 
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Table 2. Masses of neutral carbon and H 2 in high-redshift objects. 



Object 


z C o & 


M [C i] 


M^ 2 (CO) 


M^ 2 ([CI]) 


Refs. 






[M /1O 7 ] 


[M /1O 10 ] 


[M /1O 10 ] 


CO, [CI] 


IRAS F10214 


2.29 


(2.7±0.3)m- 1 


(11.3±2.5)m~ 1 


(18.8±2.4)m- 1 


1,2 


H1413+117 


2.56 


(s.i+i^m- 1 


(44.3±0.7)m- 1 




3,2 


SMM J14011+0252 


2.57 






(25.9±4.3)m- 1 


4,2 


APM 08279+5255 


3.91 


(4.4±0.6)m- 1 




(26.7±3.7)m- 1 


5,6 


PSS 2322+1944 


4.12 


(3.8±0.6)m- 1 


(12.5±5.3)m- 1 


(25.2±3.7)m- 1 


7,8 



a Redshift of CO emission. 



Note. — m is the gravitational lensing magnification factor 
References - (1) Solomon et al. 1992; (2) Weiss et al. 2005; (3) Weiss et al. 2003; (4) Downes 
& Solomon 2003; (5) Lewis et al. 2002; (6) This paper; (7) Carilli et al. 2002; (8) Pety et al. 
2005 
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Fig. 1. — Images of the [CI] 3 Pi— 3 Po and 100.2 GHz continuum emission in 
APM 08279+5255 over the full observed velocity range. The contour intervals are -2, 2, 
4, 6 and 8xer (0.36 mJy beam -1 ). The cross marks the position of the CO peak position 
from D99, which is offset (1'.'19, -0"36) from the phase center (08 h 31 m 41 8 .57, 52°45 m 17 s .7). 
The filled ellipse in the top left panel shows the 7" 2 x 6'.'0 PA 89.0° synthesized beam. The ve- 
locities indicated in the upper right corner in each panel are relative to z=3.911, determined 
from CO observations. . 
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Fig. 2. — Spectra of [CI] 3 Pi— 3 Po (top) and H 2 lio— loi (bottom) at the position of peak 
continuum intensity. Dotted lines show the fitted continuum levels, calculated from the 
intensity in the channels marked by the hatched regions. The dashed line shows a Gaussian 
fit to the [CI] 3 Pi— 3 P emission. 



